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- - ■ ABSTRACT 
\Q ■ 

O ■ 

o : 

CN . We present Keck Interferometer observations of TW Hya that spatially resolve 

its emission at 2 /im wavelength. Analyzing these data together with existing 

1—5 ! A-band veiling and near-infrared photometric measurements, we conclude that 

CN \ the inner disk consists of optically thin, sub-micron-sized dust extending from 

\ ~ 4 AU to within 0.06 AU of the central star. The inner disk edge may be 

magnetospherically truncated. Even if we account for the presence of gas in the 

fC) \ inner disk, these small dust grains have survival times against radiation blow- 

out that are orders of magnitude shorter than the age of the system, suggesting 

O ! continual replenishment through collisions of larger bodies. 

*sO . 

! Subject headings: stars: planetary systems: protoplanetary disks — stars: pre- 

main-sequence — stars: individual (TW Hya) 

O 

a: 

1. Introduction 

> ■ 
• i— i . 

TW Hya is a nearby (~ 51 pc; Mamajek 2005), young (~ 10 Myr; Webb et al. 1999) star 
surrounded by an accretion disk that evinces a large inner hole as judged from the observed 
spectral energy distribution (SED; Calvet et al. 2002). Unusually low excess emission at 
wavelengths A < 10 /im can be modeled with an optically thick disk whose inner edge is 
located ~ 4 AU from the central star. While observations of 10 /im silicate emission (Sitko 
et al. 2000; Uchida et al. 2004) together with non-zero excess at 2 /im (Johns-Krull & Valenti 
2001) suggest the presence of at least some dust grains with sizes less than a few microns 
at stellocentric distances R < 4 AU, this inner disk material appears optically thin, and has 
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been estimated to constitute less than a lunar mass (Calvet et al. 2002). Detection of warm 
gas (Herczeg et al. 2004; Rettig et al. 2004) and accretion signatures (Muzerolle et al. 2000; 
Alencar & Batalha 2002) confirm that the region inside 4 AU is not devoid of material. 

At R ~ 4 AU, dust temperatures (~ 100 K) are substantially lower than sublimation 
temperatures for silicate grains (> 1500 K; e.g., Pollack et al. 1994). This suggests the 
optically thick outer disk is truncated at 4 AU by a mechanism other than dust sublimation. 
Large holes inferred from SEDs are commonly attributed to planets, which may clear gaps 
about their orbits. A planet impedes accretion of material outside its orbit, while inner disk 
material is free to drain onto the central star (Goldreich & Tremaine 1982; Bryden et al. 
1999; Rice et al. 2003). However, a viscous outer disk causes inward migration of planets 
and their associated gaps (e.g., Lin & Papaloizou 1986; Ward 1997). Thus, inner holes would 
be filled in on the viscous timescale unless the disk is less massive than the planet, in which 
case the timescale is lengthened by the mass ratio between the planet and the disk (Chiang 
2003). Since the outer disk of TW Hya is massive (> 0.1 M ; Weinberger et al. 2002; Wilner 
et al. 2005), planets are unlikely to preserve inner clearings over the lifetime of the system 
(~ 10 Myr) unless the outer disk is unusually inviscid (a < 10~ 5 ; Shakura & Sunyaev 1973). 

An alternative explanation for SED-inferred holes is that dust grains have grown larger 
than a few microns, depleting the population of small grains that would produce the near- 
IR emission. Scattered light and long- wavelength emission from the outer (R > 10 AU) 
disk of TW Hya suggest substantial grain growth, up to cm sizes (Weinberger et al. 2002; 
Wilner et al. 2005), supporting the hypothesis of grain coagulation. The small amount of 
sub-micron-sized dust required to explain emission at A < 10 /im and the spectral shape 
of the 10 /im silicate feature (Calvet et al. 2002; Uchida et al. 2004) may represent the tail 
at small sizes of a grain size distribution that peaks at sizes much larger than a micron. 
Furthermore, as we discuss below, small dust grains are short-lived in the TW Hya disk 
and demand continual replenishment, possibly from collisions of larger parent bodies. The 
population of sub-micron-sized dust in the inner disk may not be primordial. 

Here we present observations with the Keck Interferometer that spatially resolve the 
inner disk around TW Hya for the first time. Previous observations at sub-millimeter to 
centimeter wavelengths spatially resolved emission at larger radii (> 10 AU) and enabled 
powerful constraints on the outer disk structure and dust properties (Qi et al. 2004; Wilner 
et al. 2005). With our near-IR interferometric observations, we extend this analysis to the 
inner disk. By combining spatially resolved measurements with spectral information, we 
determine the radial distribution, temperature, and approximate grain sizes of dust. We 
confirm that the inner disk is populated by small amounts of sub-micron-sized dust, and 
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show that the inner radius 1 of this optically thin disk occurs further from the star than 
previous spatially unresolved observations imply. 

2. Observations 

We observed TW Hya near i^-band (\ K = 2.14 //m; AA = 0.3 fim) with the Keck 
Interferometer (KI) on April 21, 2005 (UT 20050421). The KI is a fringe-tracking Michelson 
interferometer that combines light from the two 10-m Keck apertures and provides an angular 
resolution of < 5 mas (Colavita & Wizinowich 2003; Colavita et al. 2003). A single 130s 
observation of TW Hya was obtained between observations of three calibrators (HD 97023, 
HD 97940, and HD 99934). Figure 1 shows the uncalibrated visibilities of the target and 
calibrators. 

Data calibration is described by Eisner et al. (2005); here we summarize the procedure. 
We first determine the system visibility (i.e., the point source response of the interferometer) 
from the weighted mean of observations of unresolved calibrators (Boden et al. 1998). Source 
and calibrator data are corrected for detection biases (Colavita 1999) and integrated into 
5s blocks. The calibrated V 2 for the target source is averaged over all 5s blocks, with 
uncertainties given by the quadrature addition of the internal scatter and the uncertainty in 
the system visibility. The calibrated, normalized, squared visibility measured for TW Hya 
at r uv = 28 MA (u = 49.6 m, v = 36.7 m) is V 2 = 0.88 ± 0.05. 

3. Modeling 

We model our measured V 2 together with previous photometric measurements at A = 1- 
5 /im (Webb et al. 1999; Sitko et al. 2000), assuming photometric uncertainties of 10%. The 
circumstellar-to-stellar flux ratio is critical for modeling the circumstellar component of the 
visibility and near-IR SED (e.g., Eisner et al. 2004); we utilize a previous measurement of 
this ratio at 2 /im, r K = 0.07 ± 0.04 (Johns-Krull & Valenti 2001). Our model consists of 
the central star and an optically thin disk. We model the central star using a Kurucz stellar 
atmosphere with radius R* = R , temperature T* = 4000 K, mass M* = 0.7 M , surface 



1 In the remainder of the paper, the "inner radius" refers to the inner edge of the optically thin inner disk, 
not the inner edge of the optically thick outer disk. We take the latter to be located at R ~ 4 AU based on 
previous modeling by Calvet et al. (2002). 



-4- 



gravity logg (cm s 2 ) = 4. 5, 2 and distance d = 51 pc (Webb et al. 1999; Johns-Krull & 
Valenti 2001; Alencar & Batalha 2002; Mamajek 2005). 

The optically thin disk extends from an inner truncation radius, i? in , to an outer radius, 
-Rout- Since most of the near-IR emission is generated close to R m where the hottest dust 
resides, our results are insensitive to R ou t] for simplicity, we assume R ont = 4 AU. The mass 
surface density of the inner disk is parameterized as £ = E ( J R/AU)^ 3/2 , where R is the 
stellocentric radius and E is the surface density at R — 1 AU. 

We assume dust grains are of a single size and adopt a simple prescription for the 
frequency-dependent dust opacity, k v = kk^/^k) 13 cm 2 g _1 . For sub-micron-sized grains, 
we set (3 — 1 and kk = 10 3 . For larger grains with sizes ~ 10/im, we take (3 = and 
k,k — 10 2 . These choices are compatible with previous computations of opacities at A = 0.1- 
5 /j,m by Miyake & Nakagawa (1993). Our normalizations are ~ 10 2 times higher than theirs 
since our k v is the dust-mass opacity as opposed to the dust+gas opacity; i.e., the units of 
k v are cm 2 per g of dust. 

We compute the dust temperature under the assumption that the disk is optically thin: 

/ D \2/(4+/3) 

W*) = r.(^J • (i) 

The total flux of the inner disk is derived by dividing the disk into annuli, computing the flux 
for each annulus, and summing the annular fluxes. Similarly, model visibilities are computed 
for each annulus, and the visibility for the entire inner disk is given by the flux-weighted 
sum of the annular visibilities. The flux in an annulus of infinitesimal width dR is 

2n 

dF u (R) = —B v (T dust ) r„ R dR, (2) 
where B v is the Planck function and r v is the vertical optical depth, 



R 



-3/2 



t v {R) = kJL = y^jj ■ ( 3 ) 

The normalized visibility for an annulus extending from Ri to R2 = R\ + dR is given by the 
difference of visibilities for uniform disks having radii equal to R\ and R2: 



V(R) 



27rr uv (Rj-Rj) 



? , 2nr uv R 2 \ _ R j ( 2 ^ r uvRi 



Xd V Xd 



(4) 



2 While the computed \ogg w 4.3 for our assumed values of i?* and M*, we adopt logg = 4.5 since Kurucz 
models exist for this value. Our results are insensitive to this small difference. 
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Here r uv = 28 MA is the uv radius, A = Ax is the observing wavelength, and J\ is a first-order 
Bessel function. 

The total flux density from the disk and the central star equals 

Fu,tot = Fu,* + / dF u (R). 

J Ri n 

The veiling at 2 fxm is 

^ rRout 

r K = — dF K (R), 

^K,* jR in 

where F K is the flux density at A = \k- The squared visibility of the model at 2 /im is 

y2 = \ F K ^ + J^dF K (R)V(R) ^ 2 

FR,tot 

where V* — 1 is the visibility of the unresolved central star. 

We solve for the best-fit parameters of the model by computing the A = 1-5 /im fluxes, 
r K , and V 2 for a grid of values of i? in and S , and minimizing the x 2 -residuals between model 
and data. Uncertainties for best-fit parameters are determined from x 2 -error ellipses (e.g., 
Eisner et al. 2004). We do not include (3 or kk as free parameters and instead consider two 
model cases: (3=1 and k,k — 10 3 (intended to model sub-micron-sized grains), and (3 = 
and k k = 10 2 (representing grains sized ~ 10 /im). 

Scattered light from the disk is justifiably ignored in our modeling. At A = 1.1 and 
1.6 /im, the disk-scattered flux on angular scales of 0''4-4" is estimated to comprise 2.4% 
and 2.1% of the stellar flux, respectively (Weinberger et al. 2002). Given the blue color and 
roughly flat surface brightness profile of the scattered light within ~ 0!'8 (Weinberger et al. 
2002), we conclude that the ii'-band scattered flux within the 50 mas field of view of KI is 
< 1%. Therefore thermal emission, as traced by the K-band veiling (~ 7% of the stellar 
flux; Johns-Krull & Valenti 2001), dominates over any scattered emission. 

While we have assumed TW Hya is a single star, a low-mass stellar companion could 
contribute to the near-IR visibilities and SED. To the best of our knowledge, no stellar 
companions have been detected in previous HST imaging or radial velocity monitoring, and 
thus the presence of a luminous second star appears unlikely. Additional KI observations, 
or astrometric and/or radial velocity monitoring could test this possibility definitively. 



(5) 
(6) 
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4. Results and Discussion 

Figure 2 shows the best-fit models for (3 — and 1, together with the V 2 and SED data. 
Best-fit values for R{ n and So, and their la uncertainties are listed in Table 1. Our model 
can reproduce the KI V 2 measurement and near-IR SED of TW Hya if (3 — 1 and R in ~ 0.06 
AU. While models with (3 = can fit the SED data well (the larger quantity of SED data 
relative to the single V 2 measurement skews the fits accordingly), only (3=1 models can 
simultaneously reproduce our V 2 measurement. From these best-fit parameter values, we 
compute the temperature at the inner truncation radius T in , the dust mass Md ust , 3 and the 
2 (Urn vertical optical depth at the inner edge r^ in . 

The inner radius of the optically thin disk (0.06 AU) exceeds that inferred from previous 
modeling of spatially unresolved data (0.02 AU; Calvet et al. 2002). Our large inner radius 
leads to an inner disk temperature lower than expected for dust sublimation (Table 1), 
suggesting that an alternate truncation mechanism is necessary. One possibility is that the 
inner disk extends inward to the magnetospheric radius i? mag , where the ram pressure from 
accretion balances the stellar magnetic pressure. Although hot dust may still exist interior 
to this radius, its high infall velocity (e.g., Edwards et al. 1994) implies an optical depth 
orders of magnitude lower than that of dust outside -R mag . Assuming an accretion rate of 
5 x 10~ 10 M yr" 1 (Muzerolle et al. 2000), stellar magnetic field strength of 2.6 kG (Yang 
et al. 2005), and stellar parameters from §3, we compute R mSLg ~ 0.09 AU (e.g., Konigl 1991). 
This is comparable to our best-fit R{ n (with (3 = 1), indicating that the optically thin inner 
disk may indeed be magnetospherically truncated. 

The fact that models require (3 = 1 to fit the combined U 2 +SED data indicates that 
the inner disk contains a population of dust grains having sizes smaller than ~ 1 /im In 
fact, we can obtain slightly better fits to our data if we allow (3 > 1, as one might expect if 
very small grains (< 0.01 /im) were present (Miyake & Nakagawa 1993). 

Sub-micron-sized grains are quickly blown out of the inner disk by stellar radiation pres- 
sure. Gas friction mediates dust removal; unbound grains achieve a terminal outflow velocity 
equal to the product of the momentum stopping time and the net outward acceleration due 
to radiation pressure and gravity (e.g., Weidenschilling 1977). Estimates of the gas density 
are necessary to calculate the survival times of small dust grains. Using the measured column 
density and temperature of H 2 in the inner disk (Herczeg et al. 2004), accounting for the 



3 The dust mass depends sensitively on i? out and the assumed surface density profile. Estimates of Mdust 
are therefore highly uncertain. For comparison, under the assumption of a constant surface density, Calvet 
ct al. (2002) estimate a dust mass approximately four orders of magnitude higher than the values listed in 
Table 1. 



-7- 



possibility that the midplane may be up to 10 3 times denser than the warm surface (Glass- 
gold et al. 2004; Najita 2006), and assuming a hydrostatic disk, we estimate a midplane gas 
density at 1 AU of < 10~ 15 g cm" 3 . 4 The stopping time for micron-sized grains at R ~ 1 
AU is > 10 -2 yr, the terminal velocity is > 1 km s _1 , and the removal time is < 1 yr. 

Small dust grains in the TW Hya inner disk survive for < 1 yr, and are thus ephemeral 
over the age of the system. Because of the difficulty in transporting sub-micron-sized grains 
from the outer disk at R > 4 AU to i?i n = 0.06 AU, we argue that inner disk dust is 
continually re-generated, possibly by collisions of a swarm of larger parent bodies that also 
reside in the inner disk. 



5. Conclusions 

We observed TW Hya with the Keck Interferometer and found the 2 /im emission to 
be spatially resolved. We modeled the interferometric data together with previous measure- 
ments of the K-band veiling and near-IR fluxes, and inferred that the inner disk consists 
of optically thin dust extending from the edge of the optically thick outer disk (R ~ 4 AU; 
Calvet et al. 2002) to R m = 0.06 AU of the central star. This inner radius is larger than 
expected from dust sublimation; the truncation may be magnetospheric in origin. The near- 
IR emitting dust is composed of sub-micron- sized particles which are extremely short-lived; 
this dust may be replenished by erosive collisions of larger parent bodies in the inner disk. 
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Fig. I. — Uncalibrated, squared visibilities (V 2 ) measured by KI for TW Hya and three 
unresolved calibrator sources. The smaller V 2 value of TW Hya relative to the calibrators 
indicates that this source is angularly resolved. Scatter in the uncalibrated V 2 is due to 
a combination of instrumental and atmospheric effects, including variable performance of 
the Keck AO systems and phase jitter arising from atmospheric motions and instrumental 
vibrations. 
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Fig. 2. — Squared visibilities (left) and near-IR fluxes (right) for a model consisting of a 
young star surrounded by an optically thin disk, compared to the data (this work; Webb 
et al. 1999; Sitko et al. 2000). The flux of the central star is indicated by the solid gray line. 
Models with large dust grains (f3 = 0) can not fit the data well, in contrast to models with 
sub- micron-sized dust grains {(5 — 1). 



Table 1. Optically thin disk models 









-Rin 

(AU) 


-Rout 

(AU) 


So 
(g cm" 2 ) 


T- 

A in 

(K) 


M dust 
(g) 


T K,in 


p = 
p = 




1 


1.17 
0.98 


0.02 ±0.01 
0.06 ±0.01 


4 
4 


2.0^;? x 10~ 6 
6.3 ± 1.5 x 10~ 7 


1400 
1120 


1 x 10 22 
3 x 10 21 


0.07 
0.04 



Note. - - R m is the best-fit inner radius and S is the dust surface density 
at R — 1 AU. The outer disk radius R out is fixed for all models. The inner 
disk temperature T in , dust mass M dust , and 2 /iin vertical optical depth at the 
inner edge r^ iin , are computed for the best-fit values of R in and S . The value 
of Mdust depends on the assumed surface density profile, £ oc i?~ 3 / 2 ; different 
assumptions regarding this profile yield dust masses that can differ by orders 
of magnitude. 



